Do all liquids become strongly correlating at high pressure? 
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We present molecular dynamics simulations studying the influence of pressure on the correlation 
between the constant-volume thermal equilibrium fluctuations of virial W and potential energy U, 
focusing on liquids that are not strongly correlating at low pressure, i.e., do not have a WU corre- 
lation coefficient above 0.9. The systems studied are the two hydrogen-bonded liquids GROMOS 
methanol and TIT5P water, the ionic liquid defined by a united-atom model of the l-butyl-3-methyl- 
imidazolium nitrate and, for reference, the standard single-component Lennard- Jones liquid. The 
simulations were performed for pressures varying from GPa to 10 GPa. For all systems studied we 
find that the virial / potential energy correlation increases with increasing pressure. This suggests 
that if crystallization is avoided, all liquids become strongly correlating at sufficiently high pressure. 



The properties of strongly correlating liquids were re- 
cently discussed in several papers These liquids 
by definition exhibit strong correlations between their 
constant-volume equilibrium fluctuations of the poten- 
tial energy U and the virial 0,1 W = -1/3 £.r, • 
V r< {/(ri, rjv), where is the position of particle i. 
Recall that, if p is the pressure, V the volume, N the 
number of particles, and T the temperature, the average 
virial (W) gives the configurational contribution to the 
pressure 



pV = Nk B T + (W) 



(1) 



If A denotes the instantaneous deviations from equilib- 
rium mean values, the WU correlation is quantified by 
the correlation coefficient R defined by (where angular 
brackets denote NVT ensemble averages, i.e., averages 
at constant volume and temperature) 



R = 



(AWAU) 



V((AW) 2 )((AC/)2) 



(2) 



Perfect correlation gives R = 1 and strongly correlating 
liquids are defined |l| by R > 0.9. 

The computer simulations of Ref. [l| indicate that the 
correlation coefficient R tends to increase at increasing 
pressure, but no systematic studies have been carried out 
of the effect of pressure on the correlation. The simu- 
lations of Ref. 1 showed that van der Waals type liq- 
uids and metallic liquids are generally strongly correlat- 
ing. In contrast, liquids composed of molecules whose 
interactions have competing or directional interactions 
are generally not strongly correlating. The latter classes 
of liquids include the hydrogen-bonded liquids, the co- 
valently bonded liquids, and the strongly ionic liquids. 
Since previous works indicated that R increases at in- 
creasing pressure, one may ask whether all liquids become 
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strongly correlating at sufficiently large pressure. This is 
the question addressed in the present brief report. 

Simulations of four different model liquids were per- 
formed with NVT molecular dynamics using the Gro- 
macs package For each model samples of different 
densities were created and mixed during 1 ps (argon 
units) at a high temperature, followed by a ramping down 
to the desired isotherm. Here the systems were equili- 
brated at constant temperature during 10 - 500 ns. The 
data for each state point shown below represent an av- 
erage taken over five statistically independent samples, 
with a sampling frequency of 0.2 ps and production runs 
of length 10 ns. The following systems were studied: 
1) The single- component Lennard- Jones liquid defined by 
the pair potential Vlj{t) = 4e [(a/r) 12 — (cr/r) 6 ]. This 
system serves as a reference strongly correlating liquid. 
The results reported below refer to standard argon units 
(a = 0.34 nm, e = 0.997 kJ/mol). Samples consist- 
ing of N = 864 particles were studied. 2) Methanol: 
The GROMOS force field was used which is com- 

posed of three sites representing, respectively, the methyl 
group, the oxygen atom, and the oxygen-bonded hydro- 
gen atom (H). The masses are, respectively, 15.035 u, 
15.999 u, 1.008 u; the Coulomb interactions are given 
by the following charges: 0.176 e, -0.574 e, and 0.398 e. 
The sites interact with sites on other methanol molecules 
by additional Lennard- Jones interactions with the con- 
stants e mm = 0.9444 kJ/mol, e o = 0.8496 kJ/mol, 
€ M o = 0.9770 kJ/mol, a MM = 0.3646 nm, a 00 = 0.2955 
nm and omo — 0.3235 nm. The van der Waals interac- 
tions are cut off smoothly between 0.9 nm and 1.1 nm. 
The M-0 distance is fixed at 0.136 nm, the O-H distance 
at 0.1 nm, and the M-O-H bond angle at 108.53°. Sam- 
ples consisting of TV = 1728 molecules were studied. 3) 
TIP5P: In this water model 0] each water molecule is 
described by five sites: one site represents the oxygen 
atom (O), two sites represent the hydrogen atoms, and 
two sites locate the centers of negative charge that cor- 
respond to the oxygen lone-pair electrons. The potential 
parameters and charges used are the same as in Ref. 
Sample consisting of TV = 512 molecules were studied. 
4) [BMIMJ+ [N0 3 ]~ : A united-atom model of the ionic 



liquid l-butyl-3-methyl-imidazolium nitrate Q based on 
the GROMOS H force field. The same parameters were 
used as in Ref. @- 

As an example of WU correlations Fig. Q] shows the 
equilibrium fluctuations as a function of time of the 
T1P5P water model's normalized virial and potential en- 
ergy. Figure [TJa) gives data from a simulation at zero 
pressure at T — 475 K. The fluctuations of W and U 
are rather uncorrelated (R = 0.18). At lower tempera- 
tures the correlation is even lower; indeed near the den- 
sity maximum the correlation is close to zero [1] . Figure 
[ljb) gives data from a simulation at p = 8 GPa at the 
same temperature. Here the correlation is significantly 
larger (R = 0.64). 




FIG. 1: Time series from simulation of 512 molecules of 
TIP5P water in the NVT ensemble at two different pressures. 

(a) At zero pressure the correlation between normalized fluc- 
tuations of the virial, AW (t) / ^ {(AW) 2 } and that of the po- 
tential energy, AU(t)/y/ ((At/) 2 ), is weak, with a correlation 
coefficient of R — 0.18. As shown in Ref. Q] this low cor- 
relation is related to the existence of a density maximum at 
lower temperature where the WU correlation is almost zero. 

(b) At the pressure 8 GPa the WU correlation is considerably 
stronger (R — 0.64). 

To systematically investigate the influence of pressure 
on the WU correlation we calculated the correlation co- 
efficient as a function of pressure along isotherms for the 
four liquids. Figure [2] shows that the correlation increases 
with increasing pressure for all systems. The rather low 
correlation in the case of the single component Lennard- 
Jones (SCLJ) liquid at T=310 K reflects the fact that 
only the last three points stem from liquid-state simula- 
tions. 

In Fig. [3] the correlation coefficient was plotted in- 
stead as a function of the relative volume change, AV = 
(Vq — V)/Vq, where Vb is the highest volume at the given 
temperature corresponding to the lowest pressure of the 
simulation. In three cases the lowest pressure was around 
1 bar, i.e., effectively close to zero, but for Methanol the 
lowest pressure was of order 0.1 GPa. Water crystallizes 
upon compression before it reaches the correlation coef- 
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FIG. 2: The WU correlation coefficient R plotted as a func- 
tion of pressure along isotherms for the following systems: 1) 
The standard, single-component Lennard- Jones liquid (two 
isotherms), 2) the ionic liquid [BMIM]+ [NO3]" , 3) methanol, 
and 4) the TIP5P water model (two isotherms; the last point 
represents a crystallized sample). In all cases the correlation 
increases with increasing pressure. Data were taken from 10 
ns of simulations of each liquid. 
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FIG. 3: The correlation coefficient R plotted against the vol- 
ume decrease relative to the volume Vb at the lowest pressure 
of the given isotherm. The higher density isotherm (450 K) 
of TIP5P water shows stronger correlation than its less dense 
counterpart (475 K) at the same relative volume change. 



ficient R > 0.9 that defines a strongly correlating liquid. 

In summary, all liquids studied show increasing virial / 
potential energy correlations as pressure increases. These 
simulations indicate that if crystallization is avoided, all 
liquids become strongly correlating at sufficiently high 
pressure. 
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